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Abstract 

We examined offsets between Hn regions and molecular clouds belonging to spiral arms of a late type 
spiral galaxy NGC 4254 (M99). We used a high resolution 12 CO (J=l-0) image obtained by Nobeyama 
Millimeter Array (NMA) and an Ha image. We derived angular offsets (9) in the galactic disk, and found 
that these offsets show a linear dependence on the angular rotation velocity of gas (^g)- This linear 
relation can be expressed by an equation: 9 = (Qq — Op) -ina, where fip and t^ a are constant. Here, ftp 
corresponds to the pattern speed of spiral arms and t^ a is interpreted as the timescale between the peak 
compression of the molecular gas in spiral arms and the peak of massive star formation. We may thus 
determine fip and t^ a simultaneously by fitting a line to our 9 — VLq plot, if we assumed they are constant. 
From our plot, we obtained t Ha = (4.8 ± 1.2) x 10 6 yr and Q P = 26±g° km s" 1 kpc" 1 , which are consistent 
with previous studies. We suggest that this 9 — fio plot can be a new tool to determine the pattern speed 
and the typical timescale needed for star formations. 
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1. Introduction 

1.1. Pattern Speed 

Pattern speed (fip) is defined as an angular rotation 
velocity of spiral pattern. According to the spiral density 
wave theory (Lin & Shu 1964), flp determines the location 
and even the existence of resonances, which would greatly 
affect the morphology and kinematics of a galaxy. The 
results of numerical simulations also show the large de- 
pendence of spiral structure on f2p (Sempere et al. 1995). 
Despite its importance on the study of galaxies, it can- 
not be determined directly from observations, since the 
pattern structure is not a material but a density wave 
propagating through the disk. 

Several methods have been proposed to date for 
the determination of pattern speed. Tremaine & 
Weinberg (1984) presented the method using the con- 
tinuity equation for the surface brightness of galaxies. 
Canzian (1993) showed that the residual velocity field, ob- 
tained by subtracting the axisymmetric component from 
observed velocity field, has a difference in pattern be- 
tween the inside and outside of the corotation. Cepa & 
Bcckman (1990) pointed out that where the arm to in- 
terarm ratio of star formation efficiency has a sharp dip 
could be thought as the corotation. These methods work 
well for some galaxies, but there are still inconsistencies 
or large uncertainties. 

1.2. CO and Ha 

In this letter, we examined offsets between Ha and CO 
arms of a spiral galaxy. We suggest that we can use them 
as a new method for the determination of pattern speed 



and typical timescale for star formation. 

The Ha emission is the tracer of Hn region, which sur- 
rounds newly born massive stars, while the CO emission 
is the tracer of molecular gas. According to the galactic 
shock wave theory (Fujimoto 1968; Roberts 1969), most 
molecular clouds arc formed and then condensed by the 
compression owing to the spiral shock wave and thus be- 
long to the spiral arm. Since galactic scale star forma- 
tions take place in these molecular arms, some molecular 
gases are used to form stars, and the others are dissoci- 
ated by UV emission from massive stars born in that star 
formation. In consequence, molecular clouds disappear 
soon after their formation and their position appears to 
be fixed to the spiral pattern, though individual molecular 
clouds move at the speed of gas. Therefore, the position 
of spiral pattern and Ha arm can be measured by using 
that of a molecular arm and an ensemble of Hn regions 
respectively. In other words, offset between the arms of 
Ha and CO represents the difference between the speed 
of gas and spiral pattern, as well as the time needed for 
star formation from molecular clouds. 

This offset has been found in many spiral galaxies (eg. 
Rand & Kulkarni 1990). We can easily see it in a B-band 
image as a displacement of an arm's peak and a dust lane, 
since the light of B-band is mainly emitted by the young 
OB star and the distribution of dust lane and molecular 
gas are almost the same. However, this offset has not been 
studied quantitatively yet, since the resolution of the CO 
image was too large to examine detail structures of spiral 
arms. 
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Fig. 1. The basic idea of our method. If we observe a face-on 
spiral galaxy at f = (the left panel) , the same galaxy will be 
observed as the right panel at t = tn a . The thick solid lines 
are molecular arms at the time t of each panel. The thick 
dashed lines in the right panel (t=tH a ) show the position of 
molecular arms in the left panel (t=0). The offset distance 
between Ho and CO arm is d, expressed in equation (1). 

2. A Method for Determining the Pattern Speed 

We assume that the pattern of spiral is rigid, and 
that the gas rotates in circular orbit. Then wc define 
a timcscale tn a as an average time for Ha arm to develop 
from a galactic-shock-compressed molecular arm. If the 
physical processes of star formation are not extremely dif- 
ferent in the spiral disk, this timescale can be regarded as 
a constant parameter representing a typical value of the 
entire disk. 

Figure 1 illustrates this idea for the inside of the CR. If 
we observe a face-on spiral galaxy at t = (the left panel) , 
at t = tna the same galaxy will be observed as the right 
panel and the offset distance between the arms of Ha and 
CO, d, can be written as 
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where vq is the velocity of gas, and vp is the velocity of 
pattern. We adopt this expression, in which d becomes 
a positive value, since most of molecules exist in a rather 
central part of the disk and thus they seem to be in the 
inside of the CR. If molecular gas is more extended and 
the outside of the CR should be taken into account, Ha 
arm will be seen on the concave side of CO arm and d will 
be negative. 

Dividing both sides of equation (1) by radius R [kpc], 
we obtain 
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where £Iq = vq /R, = v p /R, and 9 is the azimuthal 
offset. This equation shows the relation between two ob- 
servables, fie and 9, and we can rewrite it as 

6>~0.58(ft G -O P ) t Ha (3) 

where 9 is in degree and <hq is in 10 7 yr. If we assume that 
iua is constant in a certain range of radius in a galaxy, 
9 becomes a linear function of fie, since we assume the 




Fig. 2. Projected image of NGC 4254: CO contours on an 
Ha image. 



rigid pattern, in other words, constant Qp. Therefore, 
by plotting 9 against Qq and fitting them with a line, 
both Qp and tpt a can be determined at the same time. In 
addition, the linearity of this plot can verify the validity 
of the assumptions we made. Wc discuss this validity in 
section 4. 

3. Data 

We observed a SAc galaxy NGC 4254 in the 12 CO ( J=l- 
0) line using Nobeyama Millimeter Array during a long- 
term project, "Virgo high resolution CO survey" (Sofue 
et al. 2003a). The spatial resolution is about 2 arcsec, 
which corresponds to ~ 160 pc at a distance from the 
Virgo cluster of 16.1 Mpc. This resolution is small enough 
to trace the spiral arms with typical width of 1 kpc. For 
more detailed information about this galaxy and the CO 
observation, see Sofue et al. (2003c). 

Koopmann ct al. (2001) presented broadband R and 
narrowband Ha images of 63 spiral galaxies in the Virgo 
cluster. The resolution of the image is 1.2 arcsec, compa- 
rable to that of our CO data. We determined the coordi- 
nate of the Ha image using the Galactic stars seen in the 
R band image, since these two images had the same field 
of view. The typical uncertainty of this coordinate fitting 
is about 1 arcsec, smaller than the resolution of the im- 
ages. We show an overlaid image of CO and Ha in Figure 
2. 

As we see in this figure, it is nearly face-on (We de- 
rived an inclination angle as 34° for central region. See 
the following paragraph for detail.) and has well-ordered 
spiral arms (classified as arm class 9 by Elmcgrcen & 
Elmegreen (1987)), so that the spiral structure can be eas- 
ily traced. This is why we selected this galaxy to examine 
offsets between Ha and CO arms. 

Position angle (P.A.) and inclination (i) arc also impor- 
tant parameters of a galaxy. We determined them apply- 
ing the task 'GAL' of AIPS to the velocity field of CO. 
This task cuts annuli out from a velocity field, and fits 




Fig. 3. Rotation curve of NGC 4254, determined by CO ve- 
locity field. 

them with a pure circular rotation velocity field in order 
to derive P. A. and i at each radius. We averaged the de- 
rived values at 5" < r < 16", corresponding to 0.4 kpc < r < 
1.2 kpc, where the values do not change largely, and then 
obtained P. A. = 72° ± 4° and i = 34° ± 5°. In Table 1, 
these P.A. and i are listed with values from other papers. 
Although only the central region was used in this work, 
the determined values arc not so different from previous 
work, in most of which the whole image of the galaxy was 
used for the determination. 

Using the obtained P.A. and i, we calculated the rota- 
tion velocity of gas by the use of 'GAL' again. The resul- 
tant rotation curve (Figure 3) is not perfectly the same 
with that derived by the iteration method (Takamiya & 
Sofue 2002; Sofue ct al. 2003b), since they used the differ- 
ent value of P.A. and i. 

These P.A. and i were also used to deproject the im- 
age into a face-on view and to transform it into a polar- 
coordinate image, a phase diagram (Figure 4). Azimuth 
is taken to be zero at the NW side of minor axis and to 
increase clockwise, in the same direction of this galaxy's 
rotation assuming the trailing arms. A range between 
the two horizontal lines in this figure (1 kpc < r < 2.6 
kpc) show the region used in the following analysis. Filled 
boxes and circles are the peaks of Ha and CO intensity, 
respectively. See the following section for detail. 

4. Application and Result 

4-1. Deriving Offsets 

We divided the phase diagram into strips of 200 pc ra- 
dial width, averaged the intensity with respect to radius 
at each azimuth in each strip, and plotted the averaged 
intensity against azimuth. We defined an offset angle 9 
as an azimuthal angular separation of intensity peaks be- 



Fig. 4. Phase diagram of NGC 4254: CO contours on Ha 
image. Thick dashed lines indicate rough positions of two 
marked spiral arms. Filled circles and boxes are the peaks of 
CO arm and Ho arm, respectively, which are defined by the 
analysis described in section 4. Labels, from A to J, are to 
identify these peaks and offsets presented in Figure 5. The 
region between the two horizontal lines (1 kpc < r < 2.6 kpc) 
indicates where we used in the analysis. 

tween the Ha and CO arms. Since arms of this galaxy are 
not totally continuous and have some breaks, we could 
not find corresponding peaks at some radii. In addition, 
we only used the range of 1 kpc < r < 2.6 kpc, since the 
spiral feature is not clear in the inner region and the rota- 
tion curve has a large dispersion in the outer region (See 
Figure 3). In Figure 4, this range is indicated by two hor- 
izontal lines. The peak of CO and Ha arm is marked by a 
filled and labeled circle and box, respectively, with a label 
from A to J, corresponding to that in Figure 5. 

The largest factor for the error in 9 is the resolution 
of CO data, which is about 2 arcsec. Thus this error is 
written as 

A(9~8.9x (JL) 1 Meg]. (4) 
Vkpc/ 

In addition, the uncertainty of the coordinate fitting of Ha 
image and the fact that the determination of offset angles 
is somewhat subjective would also be possible sources of 
error in 9. We, however, neglected these errors for the 
simplicity of calculations, so that the error in 9 could be 
larger than the value obtained from equation (4) by a 
factor of 2 at most. 

4-2. Fitting 

We plotted the derived offsets against ttc , derived from 
rotation curve (Figure 5). Errorbars in offset angle are 
A9, calculated from equation (4) at each radius. We used 
the %-square fitting method and obtained t^ a = (4.8 ± 
1.2) x 10 6 yr and Q P = 26±6° km s" 1 kpc" 1 . The dashed 
line in Figure 5 is the fitted line, and the gradient and 
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Fig. 5. The plot of offset angles against the angular rota- 
tion velocity. Labels, from A to J, are to identify the offsets, 
corresponding to the intensity peaks in Figure 4. Errorbars 
are calculated from equation (4) and the dashed line is the 
fitted line obtained by the x _sc ;U are method. The horizon- 
tal-axis-intercept and the gradient of this line corresponds to 
Qp and tn a , respectively. 

horizontal-axis-intercept of this line corresponds to the 
resultant value of tn a and Qp, respectively. 

Kranz et al. (2001) showed that the corotation radius 
for the same galaxy is about 7.5 kpc from numerical simu- 
lations, and corresponding ftp is about 20 km s -1 kpc -1 . 
They assumed the distance to the galaxy as 20 Mpc. If we 
adjust this result to the distance of 16.1 Mpc, flp should 
be about 25 km s -1 kpc -1 . The typical timescale of star 
formation has been suggested as about 10 7 yr from the 
calculation of Jeans time in molecular clouds. Our results 
are in good agreement with these previous studies. 

4- 3. Discussion 

In our method, we assumed a pure circular rotation 
and a constant time delay for star formation from molec- 
ular clouds, and we did not do any correction for extinc- 
tion in Ha data. We however know that deviations from 
these assumptions are not critically large. The streaming 
motion and the velocity dispersion would generate some 
non-circular motion, but in a spiral disk, these are usu- 
ally about 10 km s -1 (Adler & Westpfahl 1996; Combes 
& Becquaert 1997), which is very small compared to the 
circular rotation velocity of ~ 100 km s -1 . The validity 
of the latter assumption, constant <h« , can be checked by 
the linearity of 9 — fie plot (Figure 5). Since this plot 
is well fitted with a line, dependence of tn a on the envi- 
ronment should be small and included in the error of the 
derived value. 

The extinction of Ha is surely significant to luminosity 
and this might lead us to overestimate the offset values, 
since the amount of extinction becomes larger as Hn re- 
gions get closer to molecular arms and an Ha arm thus 
seems to shift to the downstream side from where it really 
is. Gonzalez & Graham (1996) showed azimuthal bright- 
ness profiles of NGC 4254 in K s -, g-, and r s -band image 
at R=75" (Figure 20 in their paper). Though the shape 



of profile is different from one to another, the positions of 
peak brightness are well consistent among them. Hence, 
we can assume that the extinction would hardly change 
the position of Ha arms and neglect the effect of extinc- 
tion on the offset values. 

5. Conclusion 

We examined the offset between the arm of Ha and CO 
of SAc galaxy NGC 4254, derived ten offset values and 
found the linear relation between the offset angle and the 
angular rotation velocity. This linearity implies that a 
rigid spiral pattern exists and that the physical processes 
of star formation can be regarded as not to extremely 
change in the disk. We emphasize that this is the first 
work to show this relationship. 

With this relation, we derived the pattern speed f2p = 
26^g° km s -1 kpc -1 , and the time delay for star forma- 
tion = (4.8 ± 1.2) x 10 6 yr. We again emphasize that 
we can derive both flp and trta simultaneously. Since the 
obtained values are consistent with previous studies, the 
way of our analysis can be a new method for the deter- 
mination of the pattern speed and the typical timescale 
needed for star formations. Thus, this method can inves- 
tigate not only the kinematics, but also the star formation 
mechanisms quantitatively. Moreover, this method can be 
applied to many spiral galaxies, since offsets between Ha 
and CO arms have been found in many spirals. 
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with their R and Ha image of NGC 4254. H. N. is finan- 
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